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We derive the Ginzburg-Landau- Wilson theory for the superconducting phase transition in two 
dimensions and in the magnetic field. Without disorder the theory describes a fluctuation induced 
first-order quantum phase transition into the Abrikosov lattice. We propose a phenomenological cri- 
terion for determining the transition field and discuss the qualitative effects of disorder. Comparison 
with recent experiments on MoGe films is discussed. 
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I. INTRODUCTION 
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In the last decade it became well appreciated that 
^3 the superconducting transition can be tuned not only by 
H changing temperature, but also by varying some other 
I parameter, like the level of disorder, magnetic field, film 
y thickness, or doping of a high-T c material, at zero temper- 
q ature (T — 0) [jj]. This is an example of a quantum phase 
U transition [Q, where the order-disorder transition is not 
1—— 'brought about by increasing the entropy of the system, 
as in the finite temperature case, but by changing the 
^ quantum mechanical ground state in a fundamental way. 
t^J- A particularly interesting example is the superconductor- 
{SJ insulator transition in two dimensions (2D) ||, where the 
04 superconducting state at T — is destroyed by increas- 
1— 1 ing the level of static disorder in the system. For a while 
Jt has been thought that the magnetic field would have a 
^— ^ similar (but not exactly the same) disordering effect on 
the superconducting ground state Jl],^||| ; recently, how- 
ever, the experiments of Mason and Kapitulnik JcJ have 
suggested that by increasing the magnetic field at T = 
I the 2D system undergoes a superconductor-metal phase 
*^ transition, and, moreover, that this transition may ac- 
Ch iually be discontinuous. The metallic state they found 
O seems to be quite unusual, in that its resistivity is anoma- 
^ lously low, and it exhibits a very large magnetoresistance. 
!> Only at a higher magnetic field the expected 2D insulat- 
k> Ing state is recovered 

. Motivated by these intriguing results, in this paper we 
03 study theoretically the T — superconducting phase 
transition in 2D and in the perpendicular magnetic 
field. We begin by demonstrating that, due to the pair- 
breaking nature of the magnetic field, in 2D (and without 
disorder) there exists a regular Ginzburg-Landau- Wilson 
(GLW) theory for the fluctuating superconducting order 
parameter even at T = 0. Long time ago, it was shown 
by Maki and Tsuzuki that the GLW theory in the 
magnetic field and in 3D is almost regular at T = 0; the 
quartic term coefficient was singular only weakly (loga- 
rithmically) as T approaches zero. We show that in 2D 
this singularity is removed by the absence of the third 
direction along which the motion would be unaffected by 
the field. The resulting GLW theory describes the quan- 



tum fluctuations of the lowest Landau level (LLL) super- 
conducting order parameter (OP), and bears close resem- 
blance to the finite-T GLW theory for a 3D superconduc- 
tor. Based on the existing understanding of the finite-T 
Abrikosov transition in 3D, we then argue that the quan- 
tum Abrikosov transition in 2D and without disorder is 
generically first-order, and propose a phenomenological 
criterion for determining the first-order superconductor- 
metal transition field H sm . In particular, we show that 
the true transition field H sm at T — is always be- 
low the mean- field critical field Hq, and that the quan- 
tum phase transition is from the superconducting into a 
zero-temperature equivalent of the vortex liquid phase: 
a metallic phase, strongly renormalized by the supercon- 
ducting fluctuations, which we propose may be related to 
the anomalous metal observed by Mason and Kapitulnik. 

In a real experiment, of course, disorder is always 
present, and may actually be quite strong, so it can be 
expected to play an important role in determining the 
nature of the transition. In principle, quenched disorder 
lifts the degeneracy of the LLL for the order parame- 
ter, and, we argue, if strong enough could bring back the 
continuous phase transition to a glassy superconducting 
phase. We provide a qualitative discussion of this phase 
transition at strong disorder within the framework of a 
self-consistent Hartree approach. The superconducting 
transition in this approach would correspond to the con- 
densation into the lowest extended Hartree eigenstate, 
that in principle needs to be computed numerically. At 
weak disorder, on the other hand, based on the existing 
renormalization group studies of the thermal Abrikosov 
transition, we expect that the superconducting transition 
at T = remains discontinuous. 

The organization of the paper is as follows. In Sec. 0, 
we derive the effective GLW action appropriate for the 
superconductor-metal quantum phase transition in 2D 
in the magnetic field and without disorder. In Sec. 



ID, 



we show how this effective action leads to a fluctuation 
induced first-order transition, and calculate the critical 
field. In Sec. IV, we discuss the qualitative effects of 
quenched disorder on the phase transition, and in Sec. 
we relate our results to experiment and other theoretical 
studies and provide a brief summary. Technical details 



are relegated to the Appendices. 



II. QUANTUM ACTION AT T = 

We begin by considering the T = action for a system 
of 2D fermions in a magnetic held (h = c=l) 



s = J2j d2rdr 



Am 



V + t 



(1) 



where "0 = tp(r,T), fj e> an d J7i denote the spin, charge, 
and mass of fermions, /i is the chemical potential, and 
H = V x A is the external magnetic field. The inte- 
gration over the imaginary time r is over the whole real 
axis (T = 0) , and we neglect the Zeeman coupling of the 
magnetic field to spin. We assume the interaction V > 
which corresponds to the s-wave attraction. The par- 
tition function is the functional integral over Grassman 
fields tp with the action S as the Boltzmann weight. Fol- 
lowing the standard Hubbard-Stratonovich decoupling 
procedure § in the Cooper channel, the action can be 
rewritten as 
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S = J d 2 rdr 

- |A(r,r)| 2 - VVA(r,T)i>\i)\ - VvA^r, r)^;Vt ( 2 ) 



in terms of the fluctuating superconducting order param- 
eter (OP) A(r, r). (The sum over the repeated spin in- 
dices is here assumed.) Since the action now became 
quadratic in the fermionic fields, in principle they could 
be integrated out so that one is left with the action in 
terms of the OP only. Furthermore, one could attempt 
to expand the resulting action in powers of the OP, 
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J d\d1 if 2 (l,2)A t (l)A(2) 



dld2d3d4 K 4 (l, 2, 3, 4)A t (l)A(2)A t (3)A(4) 

(3) 



where we denoted 1 = (ri , t\ ) for brevity. The parti- 
tion function is the functional integral over all configura- 
tions of the fluctuating OP. The quadratic and the quartic 
terms in the action can be represented diagrammatically 
as in Fig. [I]. As well known, without the magnetic field 
these diagrams would lead to singular expressions for Ki 
and K4 atT — 0, as the system would always be unstable 
towards the formation of Cooper pairs. Finite magnetic 
field, however, acts as a pair-breaker, and one may ex- 
pect it to regularize the kernels in the GLW expansion, 
similarly as disorder does for the d-wave pairing |9|. In 
the rest of this section we will show how this happens in 
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FIG. 1. Diagrammatic representation of the quadratic S2 
and the quartic S4 terms of the action in (^|). Solid lines 
correspond to single-particle Green's functions. 



2D, and obtain the GLW theory for the superconducting 
transition at T = 0. 

In the magnetic field the momentum is no longer a 
good quantum number, but the gauge invariance has to 
be maintained in the theory. This dictates that instead 
of the standard gradient expansion, the kernels in the 
GLW action can now be expanded in powers of the co- 
variant derivative V — ie*A, with e* = 2e. Furthermore, 
assuming a weak coupling V, the transition into the su- 
perconducting state can be expected at a low magnetic 
field, so to calculate the diagrams in Fig. [j] one may fol- 
low Gor'kov and use the "semi-classical" approximation 
Jll|[ll]], for the single-particle Green's function 



a(r 1 ,r 2 )=e^ ri ' r2 )G°(r 1 ,r 2 ) 



(4) 



Here v is the Matsubara frequency. The phase is deter- 
mined by 0(ri,r2) = e f** A ■ dl, integrating the gauge 
field over the straight path from r 2 to ri, and G° is the 
propagator in the absence of magnetic field. Finally, since 
the OP feels the magnetic field, one may expand it in 
terms of the Landau levels, which are the solutions of 
the single-particle Schrodinger equation in the magnetic 
field for charge 2e. Near the superconducting transition 
only the lowest Landau level (LLL) configurations of the 
OP need to be retained in the effective action fil| , since 
the higher Landau levels are not critical, and their effect 
in principle could be only to renormalize the coefficients 
in the GLW expansion for the LLL modes. Choosing the 
Landau gauge A = H(0,x,0) to represent the magnetic 
field H = Hz perpendicular to the 2D plane, the OP may 
then be expanded in terms of the LLL functions 

00 

A(r,r) = ]TC„(T)e* fc ™Vn(20; M*) = e- H ^-^\ 

— OO 

(5) 

where k = 2ir/L (L is the linear size of the system) and 
n is an arbitrary integer. 

Going to Matsubara frequency space, the quadratic 
part of the action can be written as S2 — J dw Szfu)), 
and to the lowest order in small parameter loth we find 
(see Appendix ^) 

S a (w) = [1 - .gm(2flT ff ) + V^roMKIA^r)! 2 ), (6) 



where th = £/vf is the cyclotron time with I = \ j\j2eH 
the magnetic length for charge 2e and v-p the Fermi veloc- 
ity, g = AfV/2 is the dimensionless coupling (Af = m/2-K 
is the constant density in 2D), A u (r) is the Fourier trans- 
form of A(r, r), and Q, is the usual ultraviolet cutoff. Also 



(|A w (r)| 2 ) = J dr|A w (r) 



2tt 



(7) 



Similarly, the quartic part of the action $4 
/ dhjdi0idu>2Si(uj_ ! u>±, CJ2), where at zero frequencies [|L 
(see Appendix g) , 



S i (io = uj 1 =Lu 2 = 0)=a^(\A (r)\ 4 ) (8) 
/ith the constant a = 47r[ln(l + \/2)] 2 and the average 
(|Ao(r)| 4 ) = Jdr |A_ (r)| 4 

-7 E 5(n-m + p-q)C n (0)C* m (0)C p (0)C* q (0) 
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(9) 



Note that the coefficient before (|Ao(r)| 4 ) is finite in 
a finite magnetic field. At zero field it diverges, as well- 
known. In 3D the same coefficient is oc — ln(T/7J), and 
logarithmically divergent as T — ► 0. This is a conse- 
quence of the fact that the motion in the z direction re- 
mains unaffected by the magnetic field. In 2D, however, 
the singular behavior of the Landau expansion around 
the normal state A = for the T — action is com- 
pletely cured. 

The effective GLW action for the superconductor- 
metal phase transition in 2D at T = can thus be written 
as 

Sglw[$] = J drdr [$t|9 r |$ + a|$| 2 + /3|$| 4 ] (10) 

where $ = <E>(r, t), and the partition function is Z = 
J D [$*,$] exp(— Sglw [$])■ In principle the action con- 
tains higher powers of the field and the higher order 
time derivatives, but we make the usual approximation 
in retaining only the most relevant terms. The oper- 
ator \d T \ corresponds to |o;| in Matsubara space and 
we have rescaled the field $ = {^grn ) 1 ^ 2 A to bring 
the coefficient of the \u\ term in (jlTj) to unity. Thus 
a = [g- 1 - ln(2ftT H )]/(vW) and (3 = a/{irAf). 



III. FLUCTUATION INDUCED FIRST-ORDER 
TRANSITION 

When the fluctuations of the OP are ignored (the 
mean- field approximation), the system described by the 



action (|10|) at T = undergoes a continuous phase tran- 
sition from normal to the Abrikosov superconducting 
phase as the magnetic field is decreased, at the point 
where a = 0. The mean-field critical field is thus Hq = 
2f2 2 1 (e.v\)e~ 2 l 9 , and exponentially small for a weak cou- 
pling g. One can easily verify that this satisfies the stan- 
dard relation i?o£ 2 ~ 1 where £ is the superconducting 
coherence length. This mean-field critical field is simply 
the end (T = 0) point of the standard H C 2(T) Abrikosov 
line in the H-T phase diagram of a 2D type-II super- 
conductor. However, today it is well established both 
theoretically |l4|-|l(f| and experimentally Jl^Jl^], that at 
finite T Abrikosov transition is one of those rare cases 
where fluctuations even qualitatively alter the nature of 
the transition, turning it into a first-order in the clean 
case. This may be understood as a consequence of the 
macroscopic degeneracy of the LLL manifold, which is 
a unique feature of the Abrikosov transition, and which 
greatly enhances the effect of fluctuations. Any OP con- 
figuration within the LLL can be parametrized in terms 
of its zeroes (vortices), and the Abrikosov transition is 
equivalent to a freezing transition of a complicated classi- 
cal many body system, which is then typically first order 
pTj[ . Our quantum action at T = is almost identical to 
the finite-T GLW theory for a 3D superconductor in the 
magnetic field (except for a linear instead of a quadratic 
derivative with respect to the third coordinate) , and one 
may expect that, in this case, quantum fluctuations of 
the OP may also turn the T = transition into first or- 
der. To approximately obtain the critical magnetic field 
H sm of the first-order transition, it is convenient to work 
with dimensionless quantities and first rescale the fields, 
lengths, and time as (47rl 2 Tg/3) 4 1> — > $, r/(y / 7r£) — > r, 
and t/th — * t. The action in ( |10|) then becomes 



S, 



GLW 
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drdr 



? r $t|9 T |$ + 9a |$| 
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, (11) 



where the two dimensionless couplings are g T ^ a = 
(l/r H ,a)(7r£ 2 r i f/4/3) 1 / 2 . First, consider only the OP 
configurations that have no r dependence (u> = modes) . 
For those the first term of (]ll]) vanishes, and the ther- 
modynamics of the system would depend exclusively on 
the single dimensionless coupling g a . The theory would 
then look zero-dimensional, however, this is deceiving: 
the macroscopic degeneracy of the LLL still remains. In 
fact, without the first term the partition function with 
the action (11) would be identical to the finite-T par- 
tition function that describes the vortex-liquid-to-solid 
transition, and which is known to have a weak first-order 
transition at 170, = <7m ~ —6.7. This has been established 
in detailed Monte Carlo simulations (Tj| , as well as in the 
density- functional theory (l6) . To find the transition field 
with the full r-dependence of the OP included we assume 
that the transition is driven by the same mechanism of 
growing positional correlations between vortices, the only 
difference now being that the vortices are "straight" in 
r-direction over an imaginary time scale f , instead of tr- 



We identify this time scale with the correlation "length" 
in t direction. This ansatz works remarkably well in de- 
scribing the finite-T first-order transition line in clean 
YBCO [[f9|, and on that basis we expect it to be a good 
approximation here as well. The transition field is then 
determined by the condition 



5m = g a r 2 , 



(12) 



for f > 1 in units of th- To approximately determine 
f we will use the self-consistent Hartree approximation 
to the action in (|ll|). The correlation time is then 
defined as 



(13) 



a + (|$(r)|2>/4' 



where the quantum mechanical average appearing in ( |T3] ) 
is determined self-consistently as 



<l*(r)| 2 > = ± J 



du 



SrM+5a + <|*(r)| a >/4- 



(14) 



Solving Eqs. (U)-©, one obtains the equation for the 
first-order transition field 
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where the effective coupling satisfies 
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4a|g M |Vm»i \/2eH s 

7T fcp 



(15) 



(16) 



in the Oi — > oo limit. The parameter f2j is a dimension- 
less upper cutoff we introduced to regularize the integral 
in (14). It may be chosen to yield the correct answer 
when the next order in "frequency" (i.e., ui 2 ) term is 
kept in the effective action, for example. Since f2i ap- 
pears only under a log the result is little sensitive to its 
precise value. More importantly, it is evident that g e s is 
smaller than g which in turn implies that the true tran- 
sition field H sm is below the mean-field Ho, which would 
correspond to the = case. 

In Fig. ^ we present the H-g phase diagram for the 
2D system in the magnetic field and at T = 0. The 
superconductor-metal phase transition induced by the or- 
der parameter fluctuations is first-order. The mean-field 
result represented by the dashed line is included for com- 
parison. Her we have assumed f2 = ep, f^i = 10 4 , and 
qm = —7.0. The mean- field and the true transition field 
are indistiguishably close at very small coupling, but they 
start to differ significantly at larger couplings. In princi- 
ple, one may also expect a somewhat renormalized value 
of gM from the static value (ss —7), but the difference 
does not qualitatively alter our results. 




0.20 



FIG. 2. H-g phase diagram of superconductor-to-metal 
QPT in 2D in a magnetic field. The dashed line corresponds 
to the the mean-field second-order transition. The solid line is 
determined by solving Eqs. ([[5]) and ( |l6| ) and corresponds to 
the first-order transition induced by the quantum fluctuations 
of the superconducting order parameter. 



IV. EFFECTS OF DISORDER 

So far we have completely neglected the effects of static 
disorder on the superconducting transition. Since the in- 
clusion of disorder complicates the problem significantly, 
we will discuss its effects only qualitatively. First, it 
should be possible to model disorder by including a ran- 
dom potential term V(r)|$(r)| 2 into the OP effective 
action in (fj"l|). That such a term indeed arises in the 
OP theory was demonstrated by one of us || in a re- 
lated problem of disordered c£-wave superconductor. For 
simplicity one may assume that the random potential 
is uncorrelated in space. With this term included one 
faces an interacting theory with quenched disorder (and 
in the LLL) which is notoriously difficult to analyze. So 
in what follows we assume that disorder is strong (com- 
pared to the quartic term), so that the self-consistent 
Hartree treatment of the interaction may be a reasonable 
starting point p0| . In this spirit we replace the quartic 
term in (11) as 



4> 4 



(17) 



where the average is self-consistently determined as 

l<A,(r)| 2 



(18) 



where <p n arc the LLL cigenf unctions, and e n eigenval- 
ues of the random potential, Hartree "screened" by the 
interactions 



V{v)=g a + V{v) + \mv)\ 2 ). 



(19) 



The superconducting phase transition in this approach 
would correspond to the vanishing of the lowest eigen- 



value £„ and condensation into the corresponding ex- 
tended random eigenstate [ pl| . This transition is ex- 
pected to be continuous, but to verify this scenario one 
needs to implement the self-consistent procedure numer- 
ically, which we leave for future work. The reader should 
note however, that the main effect of disorder is to lift 
the degeneracy of the LLL manifold, and thus restore the 
possibility of a continuous transition. 

In case of weak disorder the above Hartree approxi- 
mation becomes inadequate and we may only speculate 
what happens with the transition. Our prejudice, which 
we discuss more shortly, is that the transition remains 
first order for weak disorder, and turns continuous only 
at some critical disorder strength. 



V. SUMMARY AND DISCUSSION 

Our conclusion of a first order transition in the clean 
case is in agreement with the renormalization group stud- 
ies of Brezin, Nelson, and Thiaville JTJ] and more recently 
of Moore and Newman [^2) . These authors have shown 
that due to macroscopic degeneracy of the LLL the renor- 
malization group flow in dimensions less than six is al- 
ways unstable, which is usually interpreted as a sign of a 
first-order transition. Although they studied a thermal, 
and not quantum, Abrikosov transition, their methods 
and conclusions can be readily translated to our case. 
Furthermore, our conclusion is also in agreement with the 
large- N treatment of the Abrikosov transition 23 - 2^] at 
finite temperature. As for the disorder case, Moore and 
Newman J22[ also demonstrated that with weak disorder 
the renormalization group flow remains unstable, which 
would indicate that the transition is still discontinuous. 
The case of strong disorder is outside the domain of va- 
lidity of their perturbative approach, and we feel that the 
self-consistent Hartree treatment of the interaction term 
we discussed is more appropriate. 

Our theory offers a natural explanation for the hys- 
teretic behavior observed at the T = superconductor- 
metal transition by Mason and Kapitulnik. We argue 
that it is a consequence of the fluctuation induced first- 
order transition. Also, that the non-superconducting 
state appears metallic, and not insulating, may be related 
to the fact that a weak magnetic field cuts off the weak- 
localization effects, and thus relegates the localization 
effects in 2D to lower tempereatures [^7j. The observed 
anomalously small resistivity and the large magnetoresis- 
tance of the Mason and Kapitulnik metal could be related 
to the strong superconducting fluctuations near the crit- 
ical field. Additional support for this picture comes from 
the experimental observation that the resistivity away 
from the critical field H sm vanishes as R ~ (H — Hq), 
where Ho > H sm 0. This power-law follows naturally 
from our theory as follows. Away from the critical field 
(outside the critical region dominated by the OP fluctu- 
ations) the transition appears continuous, and the resis- 



tivity should vanish according to the simple (Aslamazov- 
Larkin like) scaling [E8| 



d-2 



(20) 



where £ is the superconducting correlation length. The 
GLW action (10) for the OP at the mean-field level im- 
plies d = (since the OP is confined to the LLL which 
completely eliminates the gradient terms in 2D). At the 
mean-field level £ ~ 1/(H — Hq) 1 ^ 2 , so the power count- 
ing implies 



(21) 



in the crossover region, precisely as seen in the experi- 
ment . This, of course, ceases to be valid upon entering 
the critical region, in which, as we argued, fluctuations 
eventually drive the transition first order at a lower crit- 
ical field H sm . 

Finally, we note that the GLW theory in Eq. (10) has 
the frequency dependence characteristic of a dissipative 
system. This is analogous to what was found in the the- 
ory of disordered <i-wave superconducting phase transi- 
tion at T — Q, and is related to the pair-breaking 
nature of the magnetic field. This may provide the theo- 
retical basis for the Mason and Kapitulnik interpretation 
of the superconductor- metal transition M . 

To summarize, we have studied the T = supercon- 
ducting phase transition for a 2D system in a magnetic 
field. We derived an effective Ginzburg-Landau- Wilson 
action for the fluctuating superconducting order param- 
eter which enables one to investigate the effects of quan- 
tum fluctuations and the order of the phase transition. It 
is argued that without, or with weak disorder, the quan- 
tum superconductor-metal phase transition is of first or- 
der. For the case of strong disorder we expect that the 
possibility of a continuous phase transition is restored. 
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APPENDIX A: DERIVATION OF S 2 

In this Appendix we give the explicit derivation of 
quadratic S2 in (||). If we write S 2 = f dw S 2 (w), based 
on the diagram in Fig. [I], one has at T = (apart from 
the constant term) 

V 



[ -^dr 1 dr 2 G 1/+w (ri,r 2 )G_ l ,(r l! r 2 )At(r 1 )A a; (r 2 ), (Al) 

J 27T 



where A w is given by (^J) with C n (r) — > C n (o;) in terms of 
LLL. In the semi-classical approximation, the 2D single- 
particle Green's function G is given by (||) in which 

where e p = p 2 /2m — e F and the phase 0(r l7 r 2 ) = 
eH{x\ +X2){y\ — y 2 ) / 2 . We have assumed that the field 
A = H(0,x,Q), which leads to H = Hz. After sub- 
stitution of the above into ( |Al[ ) and carrying out the 
integrations over y\ and 1/2, one obtains 



S 2 (ui) 



V 
~2 



dx\dx 2 



dp dt 

(2tt) 2 (2tt) ; 



C TO (w)5[ij, + eii(x\ + X2) — kn] 5(kn — km) x 



+ w) - e p ] [-if - (e p - v F • t; 
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(x 2 ). (A3) 



Here we have set t = q + p and e q ~ e p — vp • t with vp 
the Fermi velocity. 

Using = JVJdeJ^^ (with M = m/2n as 

the 2D density of states) and the integral (fi is the usual 
ultraviolet cutoff) 



30 ^ r de 1 

oo 2tt J_ n [i(v +w)-e] [-iv - (e - v F • t)] 



In 



_{v F ■ t) 2 + lj 2 _ 
one obtains after some calculation 

S 2 (u) = K 2 (u J ){\^{v)\ 2 ), 
where (|A w (r)| 2 ) is defined in ffl) and 



# 2 (w) = - 



.9 
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(v F -t) 2 +^ 2 



(A4) 



(A5) 



(A6) 



Here g = A/V/2. Expanding A 2 (w) = A 2 (0) + + 
But 2 + ... , we found the zero- frequency term 



K 2 (0) = - 5 ln(2fir H ) 



(A7) 



with th = £/w F and ^ = l/y2eH, and the first-order 
coefficient 



dK 2 (uj) 



duo 



(A8) 



The combination of ( A7 ) and ( A8 ) gives the result in (jg) 



APPENDIX B: DERIVATION OF S 4 

The quartic action S4 is expressed diagrammatically in 
Fig. |l|. The similar diagram has been calculated by Maki 
and Tsuzuki for 3D and finite-T. Here we need S4 



in 2D and at T = 0. One can write the quartic term as 
S4 = J dudu)idL)23i(w, ct>i, UJ2), where at zero frequencies 

(U> = OJi = U) 2 = 0) 

V 2 f dv 



5 4 (0) = — / — dridr 2 dr 3 dr4G ! ,(ri,r 2 )G_ l/ (r3,r 2 ) x 
4 J 2tt 

G u (r 3 , r 4 )G_ I/ (n, r 4 )Aj(n) A (r 2 ) Aj(r 3 ) A (r 4 ). (Bl) 

In analogy to Appendix A, after a lengthy calculation we 
find 

5 4 (0) = i 9 y« F 2 ^ ^ C*C m C;C t R(n,m,p,q), (B2) 
where for brevity, C n = C„(0), and 



i?(n, m, p, g) = 5(n — m + p — q) 
8(1,3; 2,4) 



dx\dx 2 dx 3 dx^-. 



X\ + x 3 - x 2 - x 4 1 
Here is the azimuthal angle in 2D, 



d6 1 
1 , 27r cos 2 



X 



(B3) 



0(1, 3; 2, 4) = 1, if xi,x 3 > x 2l x 4 or x\,x% < x 2 ,x 4 
0, otherwise , (B4) 



and 

T = -(x, -^) 2 + (.x 2 



2eiT 



1 / kq , 2 

+ (Xa ) 



A- 



v 2e# 

One can further solve 



2eiT V 2eH J 



(n 2 -m 2 +p 2 - q 2 ). (B5) 



1 



, 2tt tt 



5 4 (0) = - 9 y % 2 ^W i ^: £ 6(n-m + p-q)x 



where 



Z = 



4ei? 

fe 2 
16ei? 



n,m,p,q 



1 



(B6) 



n +m +p + 9 --(n + rn + p + g) 



[(n — rri) 2 + (n — p) 2 + (n — q) 2 + 
(m — p) 2 + (m - q) 2 + (p — g) 2 ] 



(B7) 



and 



f= r^^— [°°tdt Cdu [ U dve- eHt2A 
J 2tt cos 2 9 J Jo Jo 

coah(V^HtY,). (B8) 
Here (a = tan 0) 

A = i(l - ia)(l - u) 2 + ^(1 + H« 2 + J(l + ^) (B9) 



and 

£ = 



k 



2VeH 

[(1 - ia)(l - u)(n -p) + (l + ia)v{m - q)], (BIO) 

As pointed out by Maki and Tsuzuki, due to the factor 
of e~ z in (B6), there is a dominance of n = m = p = 
q when k 2 /AeH ^> 1 (i.e., in the weak-field li mit ). In 
this limit £ w 0, and thus drops out of / in (|B8|). As 
a co nseq uence, / becomes independent of n,m,p, q and 
Eq. (B6) is simplified to 



Sa(0) = - 



1 9 2 



2Afi 



-/<|A (r)| 4 



(Bll) 



where (|Ao(r)| 4 ) is given in (||). To the leading order in 
field, we found 



/ = f(H) 



47T 



[ln(l + V2)f 



(B12) 



It is interesting to note that for 3D and finite-T, the same 
term will be / cx — \n(T/H), and divergent at T — > 0. 
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